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ABSTRACT: It is argued that a unitarity-violating but weakly CPT invariant su- 
perscattering matrix exists for leading-order large-iV dilaton gravity, if and only if one 
includes in the Hilbert space planckian "thunderpop" excitations which create white holes. 
CPT apparently cannot be realized in a low-energy effective theory in which such states 
have been integrated out. Rules for computing the leading-large- N superscattering are 
described in terms of quantum field theory on a single multiply-connected spacetime ob- 
tained by sewing the future (past) horizons of the original spacetime with the past (future) 
horizons of its CPT conjugate. Some difficulties which may arise in going beyond leading 
order in 1/N are briefly discussed. 



I. INTRODUCTION 



It has been forcefully argued [hawk] that the laws of quantum mechanics, when applied 
to black holes, imply the breakdown of unitarity, at least for observers outside the black 
hole who cannot tell whether the the information has been destroyed at a singularity, 
carried away by a baby universe, or inaccessibly stored in a remnant. The S'-matrix 
must then be replaced by a superscattering matrix which in general maps pure states 
into mixed states (but conserves probability). Detailed analyses of two-dimensional toy 
models over the last year and a half[CGHS,emod,rst,gine,lowe,pira] have greatly fortified 
the arguments of [hawk] - although there are some new wrinkles [bos]. At least in some 
models information is effectively lost. Thus the role of unitarity as a cornerstone of physics 
is in serious jeopardy. 

CPT is also considered a sacred principle of physics. However, its validity is also 
seriously threatened in the context of quantum black hole processes. Indeed, it is easy to 
see [page] that unitarity violations imply a breakdown of strong CPT. Strong CPT states 
that if a given in density matrix p^ n is mapped to a given out density matrix p^ ut ] 

Pout = 0ptw (done) 

then the CPT reverse of p^ ut is mapped to the CPT reverse of p^ ut : 

Qpf n = $e- 1 p 3 ut , (gtwo) 

where ©(O" 1 ) is the CPT map from in (out) to out (in) density matrices. Multiplying 
both sides of (gtwo) by _1 and inserting it into the right hand side of (gone) one learns 
that 

0- 1 = O^O" 1 (sept) 

is the inverse of 0. However, it is easily seen [page], that if pure states go to mixed states, 
information is lost and therefore can not have an inverse. 

Put another way, the ^-matrix for black hole processes contains an intrinsic arrow of 
time: it takes pure states to mixed states but not mixed states to pure states. This is 
inconsistent with strong CPT. 

However Wald [wald] has pointed out that strong CPT violation is really a consequence 
of our CPT non-invariant description of the dynamics, rather than of a fundamental arrow 
of time. To see this consider [sork] a ^-matrix which maps any state of a spin 1/2 particle 
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to either an up or down state, both with probability one half. Clearly this involves no 
fundamental arrow of time: one could not discern whether a movie of the scattering process 
was being run forward or backward. Physical discernable violations of CPT must violate 
the weak CPT condition. This states that the probability to scatter from a given pure 
in-state pf to a pure out-state p^L is equal to the probability of scattering from 6 -1 ^ , 



Note that (wept) and (sept) are different, (wept) is obeyed in the given example of the 
spin 1/2 particle. 

Wald has speculated that the ^-matrix for black hole formation and evaporation obeys 
the weak CPT condition. An immediate consequence of this would be that the CPT 
reverse of the final state of black hole evaporation (consisting mostly of outgoing Hawking 
radiation) usually collapses to form a black hole [wald] . 

In this paper we will show that this is not the case in a two dimensional model. CPT 
reversed Hawking radiation is below threshold for black hole formation: it is reflected 
through the origin. We strongly suspect that this is also the case in four dimensions. 

However we shall see in two dimensions that weak CPT invariance can be restored — at 
least in a sector of the Hilbert space — by including the possibility of white hole forma- 
tion/evaporation*. This process has a semiclassical description which is basically the time 
reverse of black hole formation/evaporation . Black holes lead to unpredictability because 
information is lost behind the future event horizon, while white holes lead to unpredictabil- 
ity because one must integrate over initial conditions at the past event horizon. We will 
find that the integration measure is fixed by weak CPT invariance. An intriguing represen- 
tation of the ^-matrix will be given as an S'-matrix on a (in general multiply-connected) 
spacetime obtained by sewing together white and black hole horizons of the original space- 
time with its time reverse. 

* This possibility was rejected by Waldfwald] on the basis of white hole instabilities [wbld]. 
These instabilities are suppressed in a 1/N expansion, but as shall be discussed in section 
5, they may create difficulties at finite N. 





(xwz) 



Equivalently, 



T = Q~ 1 $Q 



-1 



(wept) 
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2. BRIEF REVIEW OF THE RST MODEL 



An elegant model for two-dimensional black hole evaporation was introduced by Russo, 
Susskind and Thorlacius [rst], expanding on ideas introduced in [emod]. The RST model 
differs from the original CGHS model[CGHS] by finite counterterms which are fine-tuned 
to preserve a global symmetry. This enables an exact solution of the theory at large N. Nu- 
merical analyses of the CGHS model [lowe,pira] indicate that the models are qualitatively 
similar, despite the fine-tuning. 

The classical Lagrangian for the RST model is, in conformal gauge, 



S c l = 



n 



d 2 x 



(2e^ - ^)d+d-p 



1 N 

+ e" 2 *(A V - 4d+<pd-<P) + d+fid-fi 



i=l 



(cact) 



_ 12 _ 2 , IN 
n = Jj e + ^ + Tln 



where p is the conformal factor, <fi is the dilaton and fa are N scalar matter fields. Define* 

(odef) 
(cdef) 

In the large- N limit, with x an d O held fixed, the quantum effective action is then 
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2 ' 4'" 48' 
(j) IN 



y = — e z<p + p-- In—. 

AT p 2 4 3 
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V (-d-xd+x + d+Qd-Q + \ 2 e 2 *- m ) + i £ 0+ /^-/i 
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(nact) 



When rewritten in terms of p and 0, (nact) is seen to differ from the classical action (cact) 
by the term j^d+pd-p responsible for Hawking radiation. (The effects of ghosts may be 
ignored in the large- N limit.) 

There is a residual conformal gauge invariance in (nact) . We fix this by the " Kruskal 
gauge" choice 

x = n, (gchc) 



which implies 



j. 1 ^ N 

p = <p H — in — . 
1 Y 2 12 



(rphi) 



* Our conventions differ slightly from [rst]. They are chosen so that x an d ^ are held fixed 
as N is taken to infinity. 



In Kruskal gauge the equations of motion are simply 

8+d-Q = -A 2 , 



(oeom) 



and the constraints reduce to 



din = -T$?±, (cstr) 



where 

6 N 

T ±± = ]v E d ±fi d ±fi + *±- ( tde f) 

i=l 

The functions t±(x ± ) are fixed by boundary conditions. 
The linear dilaton vacuum solution 

-\ 2 Nx + x~ 



6 = --In 

r 2 
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O = -\ 2 x + (x- + -^P+(x+)) + \m(x+) 
X z A 

1 9 
ln[-4A z x + x , 



where 



(pldv) 



4 = fe™) 

corresponds to 

O = -A 2 £ + £~ - - ln[-4A 2 x+x"] . (o/du) 
The solution corresponding to general incoming matter from Z~ is 



(gsol) 



M(x + ) = A ( X dx + x + {T^ + - t\), (mdef) 

Jo 

P + (x+)= f dx+(Tf + -t° + ). (pdef) 
JO 



'0 

and t- = By transforming back to p, <fi variables it can be seen for large M that this 
corresponds at early times to a black hole which forms and evaporates. 

However, the late-time behavior of (gsol) is unphysical. Viewed as a function of 0, O 
has a minimum at 

^ = -2 ln 48' 



There is no real value of <fi corresponding to O < O cr . O = O cr should be regarded as 
the analog of the origin of radial coordinates and the end of the spacetime, rather than 
continuing to negative radius. Reflecting boundary conditions, consistent with energy 
conservation should be imposed. RST accordingly require 



The line O = Q cr along which the boundary conditions are imposed undergoes dynamical 
motion in the x + ,x~ plane. Of course this boundary line could be moved to a fixed 
timelike coordinate line e.g. x + = x~ by a conformal transformation. However, this would 
be incompatible with Kruskal gauge and does not simplify the analysis. 

It follows from the equations of motion that the boundary curve x~(x + ) obeys 



If d+P+ is small enough, the right hand side is positive and the boundary curve is a timelike 
line. No black holes are formed: incoming matter is benignly reflected up to future null 
infinity. A similar behavior occurs in four-dimensional general relativity in that sufficiently 
weak scalar £- waves can simply pass through the origin without collapse. 

On the other hand, if <9+P+ exceeds the critical value l/4(a; + ) 2 , the boundary curve 
turns to the right (towards spatial infinity) and becomes spacelike as in the shock wave 
geometry of figure 1. It can be seen that the scalar curvature diverges along the spacelike 
segments of the boundary curve. It is not possible to implement the boundary condition 
(rbc) along these segments. Such spacelike boundary segments necessarily bound regions 
of future trapped points where <9+Q < and <9_fi < 0, which is the interior of a black 
hole. Thus these spacelike singularities resemble in every way the singularities inside four- 
dimensional black holes. We regard the singularity as a boundary of the spacetime, and 
information which flows into this spacelike curve is lost to asymptotic observers*. 

The trajectory of a spacelike segment of the boundary curve is determined, not by 
boundary conditions, but by the initial conditions on X~ . If the incoming energy is finite, 

* We are essentially defining - as opposed to deriving - the model to lose information by 
this boundary condition. Physics at the singularity is of course not controlled by our ap- 
proximations, and other boundary conditions might be possible. For example one might 
imagine a (causality violating) conveyor belt which moves the information along the sin- 
gularity and then releases it at the endpoint. It would be of great interest to find alternate 
boundary conditions which consistently implement this idea. 




(rbc) 




(ccc) 
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the boundary curve will eventually revert to a timelike trajectory. This is the "endpoint" 
at which the future apparent horizon — the boundary dividing the regions <9+fi > and 
9-1-0 < — meets the singularity, and the black hole has evaporated to zero size. After the 
endpoint the boundary conditions (rbc) are immediately imposed[rst, AS]. This is depicted 
in figure 1. 

These boundary conditions have an important consequence at the endpointf, as follows. 
Parameterizing the boundary curve as the constraints imply 

\ 2 d-x + = -d-P- + - 1 , . (hem) 

Following the boundary curve down from above the endpoint, we see that it makes a sudden 
turn to the right (away from spatial infinity) at the endpoint. This is consistent with (bem) 
only if d-P-, or equivalently T M _ is negative and infinite at the endpoint. RST indeed 
find in the exact solution that T M _ must have a negative delta-function. Since the fa are 
zero on the boundary, this implies that t- has a delta-function at the endpoint value of 
x~ . 

A general lesson, which will be important in the following, can be extracted from the 
preceding discussion. Throwing positive energy in from X - causes the boundary curve 
to lean towards spatial infinity. Above some threshold, it actually causes it to become 
spacelike. Throwing negative energy at the boundary from T~ causes it to recede. Since 
it is receding, throwing more negative energy at the boundary has less effect. It can be 
forced to make a spacelike turn to the left only by an infinitely negative energy density. 

Returning to the discussion of the endpoint, a delta function in T M _ corresponds to a 
signal which travels from the endpoint up to J + . This signal is referred to as a "thunder- 
pop". The thunderpop is an essentially planckian object, since it involves high frequencies 
and emerges from a high-curvature region. Ordinarily in a 1/N expansion wavelengths 
are restricted to be greater then 1/N, since at shorter distances energy fluctuations be- 
come large and the 1/N expansion breaks down. Thus thunderpop type states would not 
usually be included in the Hilbert space. However here we see that they generically arise 
from the evolution of long-wavelength states, and so must be considered. This is a break- 
down of decoupling. Of course as the detailed structure of the thunderpop is sensitive to 
short distance physics, the description we give here is approximate. At the very least the 
thunderpop should be smeared over a region of size 1/N. 

f RST seem to inexplicably alter their boundary condition at the endpoint. 
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Since fa = 0, yet T__ (which should be understood as an expectation value) is non-zero 
and negative along the thunderpop, it is a "highly quantum" state of / - modes. Negative 
T M _ can arise from off-diagonal terms in the expectation value [prtr]. However there are 
restrictions on how negative T^L can be: the total energy must be positive for reasonable 
states, and a region of negative T^L can not be "too far" from a region of positive T^_. 
Thus it is not obvious that the thunderpop can really be represented as a semiclassical 
quantum state. In the appendix we demonstrate that it is nevertheless possible. 

It is important to keep in mind that, although the RST model has an actual spacelike 
singularity, this will not be important in the following because the ^-matrix is insensitive to 
all physics behind the global horizons. One can certainly imagine altering the RST model 
in such a way so as to eliminate the spacelike singularity and store the information forever 
inside the black hole, or allow it to be carried away by a baby universe. The following 
discussion would apply equally to such models. 

In conclusion, the RST model embodies all the features of black hole evaporation 
anticipated by Hawking. Black holes form and evaporate in a finite time, leaving nothing 
behind. Information is lost behind a global event horizon. The only unanticipated feature 
is the generic occurrence of a thunderpop at the evaporation endpoint. In section 4 we 
shall see this plays a key role in the restoration of CPT invariance. 



3. THE RST ^-MATRIX 

Given the quantum state on X + , it is not in general possible to reconstruct the quantum 
state on X~ . Additional information, corresponding to the state swallowed by the black 
hole, is required. One therefore expects that only a ^-matrix, rather than an S'-matrix, 
exists for the RST model. It is perhaps possible to obtain analytic formulae for the matrix 
elements to leading order in large- N using technology developed in [ginejf . However, such 
explicit formulae will not be required for our considerations, and we will content ourselves 
in this section with a general description of the large- N ^-matrix. 

Incoming states are created by acting with /-oscillators on the /-vacuum. Let us first 
consider states obtained by acting on the vacuum with order N® creation operators. It 

f This has been partially worked out for shock wave geometries by L. Thorlacius (private 
communication) . 
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follows from (nact) or (tdef) that there is no back-reaction on the metric or dilaton to 
leading order in 1/N. This is simply because "Newton's constant" is taken to zero as 1/N. 
Thus the scattering is simply obtained by reflecting the /-particles off of the boundary at 
fi = fi cr , (according to (rbc)), which in the vacuum is the timelike line 4\ 2 x + x~ = — 1. 
Clearly a unitary S'-matrix exists in this sector of the Hilbert space. 

In order to investigate black hole dynamics, one must consider states with incoming 
energies (i.e., (t¥±_)) of order one. Semiclassical configurations corresponding to collapsing 
matter are typically described by specifying c-number initial data p c (x + ), (p c (x + ), f?(x + ) 
and t+(x + ) on X~ . This description is redundant because of conformal gauge invariance. 
The semiclassical quantum state corresponding to this initial data has the simplest de- 
scription in the gauge £+ = (the gauge transformation properties of £+ are described in 
the appendix). In this gauge it can be represented on X~ as a coherent state of the form: 

\ff) = A : e -^Etii 'dx+d + f?{x+)Mx+) . | Q)) {chst) 

where the normal ordering and the vacuum are defined with respect to t+ = coordinates 
and A is a normalization factor. Its time evolution is then given - to leading order in N - 
by classical evolution of ff. In general the coherent states are an overcomplete basis, but 
in the semiclassical large- A?" limit A becomes small and they become orthogonal. 

In t+ = coordinates the Kruskal condition (gchc) will of course not hold in general. 
Rather one has 

tt = X + u, (tgc) 

where a; is a c-number solution of the free wave equation. In this t+ = gauge the operator 
constraints become 

6 N 

—d\uo — d+uid+u) + + d\Vl = — — ^ d+fid+fo. (tcs) 

i=\ 

Q is then found in terms of the physical degrees of freedom fa by integrating this equation, 
with some integration constants fixed by (oeom)*. 

If the c-number field configurations ff in (chst) are such that 

fi N 
i=l 

* This is the analog of light-cone gauge in string theory. A DDF-like description, in which 
X and fi are independent operators, is also possible[verl,deal]. 
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is order iV°, there is a leading order back reaction on the geometry, and p and are 
non-zero. If p is non-zero, then the representation (chst) of the incoming state will not 
correspond to the natural Fock basis of particles which might be seen by asymptotic inertial 
observers. The Bogoliubov transformation associated to the coordinate transformation 
from the Kruskal gauge (gchc) to p = gauge can be used to reexpress the state \ff) in 
this basis. It will not in general be a simple coherent state in this asymptotic Fock basis. 

If (TjK) is below threshold for black hole production, the semiclassical geometry is 
obtained by solving the equations (oeom) with initial conditions determined by the ex- 
pectation value of (tcs) and boundary conditions of the form (rbc) suitably modified to 
conform to the gauge (tgc). The outgoing coherent state on I + into which |/?) scatters is 
found by reflection off of this dynamically determined boundary. This will naturally lead 
to an expression in terms of oscillators which are positive frequency with respect to £+ = 
coordinates on J + . However, in the process of solving for the semiclassical geometry, the 
conformal factor p is altered on J + . Thus a Bogoliubov transformation is required to 
obtain the outgoing state in a natural Fock basis on J + . 

Let us now consider further acting on the state in (chst) with order f creation 
operators. To leading order in 1/N, this will not affect the stress tensor (texp) and the 
semiclassical geometry and boundary curve will be unaltered*. The outgoing state is then 
determined by boundary reflection and Bogoliubov transformation. In this manner the S- 
matrix for states obtained by acting with order iV° oscillators on (chst) can be constructed. 

If (T^ + ) exceeds the critical threshold, then the boundary curve will have a singular 
spacelike segment. However, the incoming state on X~ will still scatter (given the endpoint 
prescription of section 2) to a definite out-state on J + LiH + , where H + is the future event 
horizon. States obtained by acting on the incoming coherent state with order iV° f- 
oscillators will scatter to a definite pure out-state on 1 + UH + . These out-states are of the 
general form 

\Vout) = ^2^xa\x)\a) (sdec) 

x,a 

where |x)(|a)) is a state on 1 + (H + )- The "S'-matrix" for such states has the form 

(x|(ct|S , |a) = S X a a - (smat) 

* Energy conservation should be restored by including the back reaction at the next order 
in 1/N. 
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In this expression, \a) denotes an incoming state which collapses to form a black hole, 
| a) denotes a state on 7i + , and \x) denotes a state on 1 + which results from black hole 
evaporation — including a thunderpop. 

The RST dollar matrix can now be constructed by tracing over states on H + 

(0rst) x y b a = S^ a b S xa a . (drst) 

first conserves probability 

($rst) x X b a = h a - (Peon) 

However, it is not the product of two unitary matrices, and pure states will evolve into 
mixed ones. 



4. CPT INVARIANCE 



It is easy to see that $ rs t is not, on its own, weakly CPT invariant. The CPT transform 
of the incoming states which collapse to form a single black hole will never appear as an 
out-state on J+: 

(x|0|a) = 0. (hrth) 

This is because the out-states all contain planckian thunderpops where (T^L) diverges. 
Also (T+ + ) always exceeds the critical value somewhere on J~ (otherwise a black hole 
would not form), but (T^L) will generically not exceed the critical value on J + (otherwise 
a white hole would be in the past of X + ). The orthogonality condition (hrth) between in 
and out Hilbert spaces immediately precludes the possibility of weak CPT invariance[wald]. 

Stated in this manner, the remedy is obvious. One must include in the incoming 
Hilbert space thunderpop states which create white holes, as illustrated in figure 2. One 
may then construct the large- N ^-matrix 

ficpt = first + ®$rst®- (sccpt) 

The first term acts on states which collapse to form black holes, while the second acts on 
thunderpop states which create white holes. It is easy to check, in bases for which O -1 = 
O r , that $ cp t obeys the weak CPT condition (wept) as well as probability conservation 
(peon). 
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Physically, the second term in (sccpt) can be viewed as follows. A thunderpop carrying 
negative, divergent (7++) arrives at the boundary curve, forcing it to turn left (away 
from spatial infinity) and become spacelike. The future of this point will then depend on 
boundary conditions along the null surface H~ . One must sum over all possible states on 
H~ , and construct a ^-matrix by tracing over these states. The measure for this trace is 
fixed by weak CPT invariance. 

Note that the existence of the thunderpop appears crucial to the restoration of CPT 
invariance. If the black holes just quietly disappeared into the vacuum, it is hard to see 
how the CPT reverse of the outgoing state could form a white hole. The thunderpop 
is required to nucleate the white hole into which the CPT-reversed Hawking radiation 
subsequently falls. This makes it difficult to see how CPT could be realized in a long- 
wavelength effective theory in which thunderpops are integrated out. It would certainly be 
interesting to understand the analog of thunderpops in four dimensions. CPT would seem 
to require that there is something which can nucleate a white hole. Perhaps planckian 
negative energy densities are necessary for four as well as two-dimensional white hole 
nucleation. 

The ^-matrix is naturally represented in the language of functional integrals. The 
scattering amplitude from a coherent state |/(J - )) on J~ which collapses to form a black 
hole to an out-state \f"(H + ) : f'(l + )) is given by 

rf(l + )J"(H+) 

S(f(H + )j\l+yj(l-))= / Vfe~ iS , (sfnt) 

where the boundary conditions are those appropriate to coherent states. The primes on the 
boundary data / indicate that they are generally different functions. The matrix element 
$ rs t between an initial coherent density matrix |/)(/| and a final coherent density matrix 
|/') (/'| is then obtained by multiplying (sfnt) by its complex conjugate, and functionally 
integrating over f"(H + ) (i.e., tracing over states on the horizon) 

Vfe~ lS / Vfe lS . 

The integration over common values of boundary data /" on H + and H + sews together 
the two path integrals along the horizons. This may be represented as a single functional 
integral on the unusual semiclassical geometry depicted in figure 3 (corresponding to white 
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and black hole spacetimes sewn together at the horizons): 

r f(I+)J'{i-) 

$[~f(i + ), /(J-)] = / ' A Vfe~ iS . (ggg) 

Jf(l-)J(1+) 

It should be understood in (ggg) that the time flows forward from X + to X~ and the sign 
of the action is reversed in the double of the original spacetime . 

Let us summarize the rule for computing the large- N dollar matrix element $[f'(I + ), /'(X + ); /(X - ), / 
for the case in which |/) and |/) both collapse to form a black hole, and \f') and \f') are 
both outgoing thunderpop states. The first step is to determine the semiclassical geometry. 
To do so, construct the stress tensor corresponding to /(X - ), and determine the incoming 
geometry by solving the constraints. Then solve the semiclassical equations of motion with 
this initial data. This will determine the full semiclassical spacetime geometry and field 
configurations, including f"(H + ) and f'(T + ). The outgoing state is (in t+ = gauge) 
the coherent state associated to this semiclassical data. The large- A?" approximation to the 
^-matrix vanishes unless the semiclassical field configuration f'(l + ) so obtained agrees 
with its argument because coherent states constructed with different c-number functions 
are orthogonal at large N. Next, solve the constraints using /'(X + ) to get (CPT reversed) 
initial data for X + . This initial data will produce a white hole. Use the CPT reverse of 
the final data f"(H + ) obtained previously on H + as initial data on the white hole horizon 
on H + . One may then compute /(X - ) as the CPT reverse of the final data in this space- 
time. The semiclassical coherent state ^-matrix element vanishes if the semiclassical data 
thereby obtained does not agree with its argument. 

^-matrix elements of states obtained by acting on these coherent states with order 
iV° oscillators are now simply .S-matrix elements (remembering to perform appropriate 
Bogoliubov transformations) in the sewn spacetime of figure 3. 

So far we have discussed processes involving one black hole or one white hole. The 
construction of the ^-matrix for two black holes or the CPT reverse, two white holes 
proceeds similarly. Now there are two spacelike segments along which the two spacetimes 
must be sewn together. 

Clustering however requires that we include initial states with both black holes and 
white holes. This can result, after sewing, in a spacetime with closed timelike loops as 
in figure 4. New features, which we have not yet analyzed, may arise in computing the 
^-matrix from S'-matrix elements on such spacetimes. 
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This prescription for computing the ^-matrix formally generalizes to four dimensions: 
Form the double of a spacetime by sewing its black (white) hole horizons to the white 
(black) hole horizons of its CPT conjugate. The ^-matrix is then obtained from the S- 
matrix on the doubled spacetime. This picture suggest a natural way of including the effects 
of virtual black holes: they are simply wormholes which connect the original spacetime to 
its double. For the usual reasons [worm], this will not lead to quantum incoherence, but 
rather will turn coupling constants into dynamical variables. The only slight difference 
is that the wormhole a-parameters label ^-matrix rather than S'-matrix superselection 
sectors. This leads us to question the assertion [hwktwo, bps] that information loss from 
real black holes inevitably implies that virtual black holes lead to tiny unitarity violations 
in a long distance effective theory. 



5. SOME PROBLEMS 



In this paper we have described the black and white hole superscattering matrix to 
leading order in a 1/N expansion. Going beyond leading order may not be straightforward 
for several reasons. 

The first is that the RST boundary condition - which entered crucially into our discus- 
sion - has only been fully defined at the semiclassical level. Heuristically one might hope 
to define the full quantum theory by restricting the range of functional integration of the 
field O. However it is not clear what this means in practice, or if the resulting theory will 
conserve energy. Furthermore virtual fluctuations of Q about its vacuum value will lead 
to internal boundaries, or holes in spacetimes. These difficulties will certainly arise at the 
non-perturbative level, and perhaps even in 1/N perturbation theory. 

A second problem is the white hole instabilities discussed in [wbld]. This problem has 
also been stressed - in a slightly different form - in more recent papers [brt,verl]. Consider a 
particular white hole formation/evaporation process obtained as the time reverse of black 
hole formation/evaporation, projected on to some particular pure state on 7i~ and Z~ . 
Now consider adding a single incoming / quanta to the state on X~ somewhere near or 
shortly after the incoming thunderpop. It is straightforward to show that the energy {i.e., 
of this state is blueshifted by an amount of order e M l^, where M is the white hole 
mass. The energy of the quanta is then of order e M ^/N. 

In the large- N limit the back reaction of this quanta can be neglected. However if N 
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is large but finite and M is greater than IniV, the back reaction is not small. The 1/N 
expansion presumably breaks down for such matrix elements at such large values of M. 

Physically one expects, for sufficiently large M, that the blueshifted quanta will lead 
to a black hole after the white hole. Thus beyond the 1/N expansion one does not expect 
in general to always be able to determine the numbers of black and white holes from the 
state on X - alone: knowledge of the state on 7i~ is also required. This will bring new 
features into the analysis, and we do not know how one develops a systematic perturbation 
expansion in this situation. 

In conclusion, the inclusion of white holes restores CPT to leading order in a 1/N 
expansion for processes with one white hole or one black hole. The general case remains 
unresolved. We hope to have convinced the reader that it is an interesting problem. 
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APPENDIX 



In this appendix we will describe the RST thunderpop [rst] as a quantum state of 
/-particles. 

For an incoming shock wave \(T+ + — t+) = M5(x + — 1) , the field Q in Kruskal gauge 
is given by 

fl = -\ z x + (x +-^) + - 

1 A A (gsoip) 

- -ln[-4A 2 x + x _ ], 

in the region above the shock wave at x + = 1, where O > O cr and prior to the thunderpop. 
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The thunderpop orignates at the black hole endpoint 

- A = f - M A 

Xs ' Xs J V A3(l-e- 4M /A)'4M 



(e 4M / A -l)), 



(xs) 



where the apparent horizon and singularity meet. For x > x s , the solution is joined onto 
the shifted linear dilaton vacuum 



o , Ml o , M 

O = -\ 2 x+(x- + jL)-- ln[-4X 2 x + (x- + ^ 



(gip) 



O is continuous along the thunderpop at x = x s . (glp) and (gsolp) obey the constraints 
(cstr) with 

1 



t- = 



_2 



Ax' 

t = -^4t^-(1 - e-± M ' x ) 2 5{x- - x-) 



t- = 



AM 
1 



A(x- + §)1 



x x q •, 



x x s , 



X > x s . 



(zzz) 



The metric corresponding to (glp) and (gsolp) is asymptotically flat on J + . Quantum 
states on 1+ are best described in terms of creation and annihilation operators defined 
with respect to the natural flat coordinates. Defining 

\ 3 x~ +M~ 



V = - In 
y + = ln(\x+), 
one finds that asymptotically on X + 



X 3 x7 + M 



(xxyy) 



ds z 



-dy dy~ 



(syy) 



In these coordinates, the radiation flux on T + is given by To determine t_, note that 
the constraints in a general gauge are 



N 

2(e- 2 ^ + =\fl 9-fid-fi + ^(t-~ d-pd-p + dip), 



(gcst) 



i=l 



where V- is the covariant derivative. Both sides of this equation must transform covari- 
antly under a conformal coordinate transformation x~ = x~(x~). This requires that 



{d-x-ft- = t-- {d-x-) l / 2 dl(d-x-)- 1 / 2 . 
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(ttr) 



One then finds that in y coordinates 

_ 2eMe 4M/A - 1) + e 2 Me 4M / A - l) 2 
4(l + e^(e 4M A_i))2 

t- = ~ e- 4M / A )%-) 



y~ <o, 



o 



2/" 



= 0, 
> 0. 



(2/2/2/) 



A thunderpop state \T) clearly can not be described as a coherent state in the Fock 
basis associated to y coordinates because the energy density of such states is non-zero only 
when (/j) is non-zero. To represent \T) as a coherent state, we must find a coordinate 
system in which t- vanishes. This is accomplished by solving (ttr) with the left hand side 
set to zero and £_ given by (yyy). One finds that t- vanishes in z coordinates 



= 2/ + , 
= - In 



z = 



e -y--4M/\ +1 _ e -4M/X 
4y - e -4M/A 



4 + (1 - e~ 4M / x )y- 
A thunderpop may then be represented as the z vacuum 



2T<0, 



y~ >0. 



(zzyy) 



\T) = |0 



Z 1 



(tpop) 



where \0 Z ) is annihilated by modes of the / field which are negative frequency with respect 
to z~ . To represent this state in the y Fock basis - which corresponds to particle states de- 
tected by inertial asymptotic observers - one must perform the Bogoliubov transformation 
corresponding to (zzyy) . A nice discussion of these transformations in two dimensions can 
be found in [gine]. 

Note that an isolated thunderpop - unaccompanied by nearby positive energy density 
- can not be represented by a quantum state as in (tpop) . The analog of the z coordinates 
in which t- vanishes would not cover all of X + . Thus the global existence of z coordinates 
is an important consistency check for the RST model. 

Finally, we would like to remind the reader that the state (tpop) is not the full outcome 
of shock wave collapse. The state on J + will be a mixed state, or density matrix. To 
find that density matrix one must first compute the Bogoliubov transformation from the 
incoming pure state on Z~ to the outgoing pure state on H + U X + , and then trace over 
states on H + . (tpop) is presumably one of the pure states which appears on J + with some 
probability. 
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FIGURE 1 . A shock wave collapses to form a black hole. The final 
mixed state consists of outgoing Hawking radiation and a thunderpop. 
The region of future trapped points is shaded. 




FIGURE 2. The CPT reverse of black hole formation/evaporation 
depicted in figure 1 is a white hole which is nucleated by an incoming 
thunderpop and subsequently evaporates. 




FIGURE 3. A doubled spacetime is formed by sewing the future horizon of figure 1 to the 
past horizon of figure 2. Superscattering matrix elements are obtained from S-matrix elements 
on this doubled spacetime. The fat arrows denote the direction of time flow. 




FIGURE 4. The doubled spacetime appropriate to superscattering processes with one 
white and one black hole. The two horizons connected by the fat arrow are identified. 



